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WEST Search History 



DATE: Monday, September 19, 2005 

... . „ Set ^ Hit 
H-de? ^ Quer* ^ 

DB=PGPB f USPT f USOC,EPAB,JPAB,DWPI, TDBD; PLUR=YES; OP=ADJ 



□ L73 



L72 and ((horizontal$3 or vertical$3 or transvers$5 or longitud$4 or parallel or ^ 
perpendicular$2 or longitudinal$2) adj4 (resist$5) with (ratio)) 

L71 and ((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel or 

□ L72 perpendicular$2 or longitudinal$2) adj4 (permeability or permeabl$3) with 2 

(ratio)) 

□ L71 5463549 15 

L57 and ((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel or 

□ L70 perpendicular$2 or longitudinal$2) adj4 (permeability or permeabl$3) with 4 

(ratio)) 

□ L69 L68 and (laminat$4) 0 

□ L68 L67 and (estima$6 or approximat$6) 1 

□ L67 L66 and (coarse or fine or water or M h20" or "h.sub.20") 1 

□ L66 L63 and (Waxman or Smits or Thomas or Stieber) 1 

□ L65 L63 and (bulk or content) 0 

□ L64 L63 and ((magnetic adj resonance) or MRI or NMR) 0 

L53 and (((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel 
or perpendicular$2 or longitudinal$2) adj 3 (permeability or permeabl$3)) with ^ 
((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel or 
perpendicular$2 or longitudinal$2) adj3 (resistiv$5))) 

L55 and (((horizontal $3 orvertical$3 or transvers$5 or longitud$4 or parallel 
or perpendicular$2 or longitudinal$2) adj 3 (permeability or permeabl$3)) with ^ 
((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel or 
perpendicular$2 or longitudinal$2) adj3 (resistiv$5))) 

L59 and (((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel 
or perpendicular$2 or longitudinal$2) adj3 (permeability or permeabl$3)) with ^ 
((horizontal$3 or vertical$3 or transvers$5 or longitud$4 or parallel or 
perpendicular$2 or longitudinal$2) adj3 (resistiv$5))) 

L59 and (((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel 
or perpendicular$2 or longitudinal$2) adj2 (permeability or permeabl$3)) with ^ 
((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel or 
perpendicular$2 or longitudinal$2) adj (resistiv$5))) 

L58 and ((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel or 
perpendicular$2 or longitudinal$2) with (permeability or permeabl$3) with 
L59 (formation or wellbore or "well bore" or well-bore or borehole or bore-hole or 5 
"bore hole" or earth or petrophysical$3 or petro-physical$3 or "petro 
physical $3" or sand or clay or shale or grain)) 



□ L63 



L61 



L60 
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□ L58 L57 and (logging) 5 

L55 and ((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel or 

□ L57 perpendicular$2 or longitudinal$2) with (permeability or permeabl$3) with 8 

(ratio)) 

L55 and ((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel or 

□ L56 perpendicular$2 or longitudinal$2) with (resistiv$5) with (permeability or 4 

permeabl$3) with (ratio)) 

□ L55 L54 and (ratio) 36 

L53 and ((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel or 

□ L54 perpendicular$2 or longitudinal$2) with (resistiv$5) with (permeability or 61 

permeabl$3)) 

n ^. - L52 and ((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel or ^ ^ 
perpendicular$2 or longitudinal$2) with (resistiv$5)) 

□ L52 L51 and (resistiv$5) 8594 

□ L51 L50 and (permeability or permeabl$3) 77352 

L49 and (formation or wellbore or "well bore n or well-bore or borehole or 

□ L50 bore-hole or "bore hole" or earth or petrophysical$3 or petro-physical$3 or 1205345 

"petro physical$3" or sand or clay or shale or grain) 

□ L49 (h°ri zonta l$ 3 or vertical $3 or transvers$5 or longitud$7 or parallel or 75 10108 

orthogonal$3 or perpendicular$3) 

LI and (formation or wellbore or "well bore" or well-bore or borehole or bore- 

□ L48 hole or "bore hole" or earth or petrophysical$3 or petro-physical$3 or "petro 683458 

physical$3" or sand or clay or shale or grain) 

□ L47 113 and (permeability or permeabl$3) 6 

□ L46 5656930 and (permeability or permeabl$3) 15 

□ L45 US20040140801A1 1 

□ L44 L43 and (permeability or permeabl$3) 5 

□ L43 3479581 11 

L37 and ((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel or 

□ L42 perpendicular$2 or longitudinal $2) with (resistiv$5) with (permeability or 5 

permeabl$3) with (ratio)) 

L40 and ((horizontal$3 or vertical$3 or transvers$5 or longitud$4 or parallel or 

□ L41 perpendicular$2 or longitudinal $2) with (resistiv$5) with (permeability or 4 

permeabl$3) with (ratio)) 

L39 and ((horizontal$3 or vertical$3 or transvers$5 or longitud$4 or parallel or 

□ L40 perpendicular$2 or longitudinal$2) with (resistiv$5) with (permeability or 36 

permeabl$3)) 

□ • L39 L38 and (ratio) 126 

L37 and (formation or wellbore or "well bore" or well-bore or borehole or 

□ L38 bore-hole or "bore hole" or earth or petrophysical$3 or petro-physical$3 or 199 

"petro physical$3" or sand or clay or shale or grain) 
— L36 and ((horizontal $3 or vertical$3 or transvers$5 or longitud$4 or parallel or 36Q 

perpendicular$2 or longitudinal$2) with (resistiv$5)) 

□ L36 L35 and ((horizontal$3 or vertical$3 or transvers$5 or longitud$4 or parallel or 1674 
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□ 
□ 

□ 
D 

□ 



□ 
□ 

□ 
□ 
□ 
□ 
□ 
□ 
□ 
□ 
D 



□ 
□ 



□ 
□ 



□ 
□ 
□ 
□ 
□ 



L35 
L34 

L33 

L32 

DB= 
L31 

DB= 



□ L30 



L29 
L28 
L27 
L26 
L25 
L24 
L23 
L22 
L21 
L20 
L19 

L18 

L17 
L16 



□ L15 



L14 
L13 



□ L12 



Lll 
L10 
L9 
L8 
L7 
L6 



perpendicular$2 or longitudinal $2) with (permeability or permeabl$3)) 
L34 and (permeability or permeabl$3) 
L33 and (resistiv$5) 

(horizontal$3 or vertical$3 or transvers$5 or longitud$4 or parallel or 
perpendicular$2 or longitudinal $2) 

(horizontal$3 or vertical$3 or transvers$5 or longitud$7 or parallel or 
perpendicular$2) 

PGPB,USPT,EPAB,JPAB,DWPI,TDBD; PLUR=YES; OP=ADJ 

20020101235 

USPT,PGPB,JPAB,EPAB,DWPI,TDBD; PLUR=YES; OP=ADJ 

6005389 and ((horizontal $3 or vertical$3 or transvers$5 or longitud$7 or 
parallel or orthogonal$3 or perpendicular$3) with (resistiv$5)) 

L28 not L26 

L27 and (Waxman or Smits or Thomas or Stieber) 
L22 and ((magnetic adj resonance) or MRI or NMR) 
L25 and ((magnetic adj resonance) or MRI or NMR) 
L24 and (coarse or fine) 

L23 and (Waxman or Smits or Thomas or Stieber) 

L22 and (bulk or content) 

L21 and (estima$6 or approximat$6) 

L20 and (induct$5) 

LI 9 and (model$4 or simulat$6) 

LI 8 and (density or porosity or permeability or bvi or irreducible or bound) 

LI 7 and (formation or wellbore or "well bore" or well-bore or borehole or 
bore-hole or "bore hole" or earth or petrophysical$3 or petro-physical$3 or 
"petro physical$3" or sand or clay or shale or grain) 

L16 and (water or fluid$4 or liquid$4 or "h20" or oil) 

((horizontal $3 or vertical$3 or transyers$5 or longitud$7 or parallel or 
orthogonal$3 or perpendicular$3) with (resistiv$5)) 

L14 and ((horizontal $3 or vertical$3 or transvers$5 or longitud$7) with 
(resistiv$5)) 

6255819 

L12 and ((magnetic adj resonance) or MRI or NMR) 

LI 1 and ((horizontal$3 or vertical$3 or transvers$5 or longitud$7) with 
(resistiv$5)) 

L10 and (bulk or content) 

L9 and (Waxman or Smits or Thomas or Stieber) 

L8 and (model$4 or simulat$6) 

L7 and (coarse or fine) 

L6 and (density or porosity or permeability or bvi or irreducible or bound) 
L5 and (estima$6 or approximat$6) 



14130 
166874 

7399585 
7403510 



3 

13 
60 
10 
21 
32 
223 
328 
359 
792 
1954 

3125 

5881 
17711 



14 
10 

12 

184 
196 
1390 
2283 
5086 
6081 
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□ 




□ 


L5 


1— 1 


T A 


□ 


L3 


□ 


L2 


□ 


LI 



L4 and (induct$5) 7901 

L3 and (formation or wellbore or "well bore" or well-bore or borehole or bore- 
hole or "bore hole" or earth or petrophysical$3 or petro-physical$3 or "petro 30322 
physical$3" or sand or clay or shale or grain) 

L2 and (water or fluid$4 or liquid$4 or "h20" or oil) 50898 

LI and (resistiv$5) 94731 

(horizontal$3 or vertical$3 or transvers$5 or longitud$7) 4225033 



END OF SEARCH HISTORY 
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Hit List 



• Clear j GenerateCollection I Print I Fwd Refs | Bkwd Refs 

llpperate OACS : ; J 



Search Results - Record(s) 1 through .10 of 10 returned. 



□ 1. Document©: US 200501 14029 Al 

Using default format because multiple data bases are involved. 

L13: Entry 1 of 10 File: PGPB 



May 26, 2005 



PGPUB-DOCUMENT-NUMBER: 20050114029 
PGPUB- FILING- TYPE : new 

DOCUMENT-IDENTIFIER: US 20050114029 Al 

TITLE: Determining water saturation for oil bearing thin-bedded formation having 
anisotropic resistivity 

PUBLICATION-DATE: May 26, 2005 



INVENTOR- INFORMATION : 
NAME 

Clavaud, Jean-Baptiste Nicolas 
Boyd, Austin Joseph 



CITY 
Houston 



STATE COUNTRY RULE- 4 7 
TX US 



Ridgef ield CT 



US 



US-CL-CURRENT: 702/7 



IHBIEll 



!;:iClassfti$a'Mf& 



□ 2. Document ID: US 20040140801 Al 

L13: Entry 2 of 10 File: PGPB 



Jul 22, 2004 



PGPUB-DOCUMENT-NUMBER : 20040140801 
PGPUB-FILING-TYPE: new 

DOCUMENT-IDENTIFIER: US 20040140801 Al 

TITLE: Combined characterization and inversion of reservoir parameters from 
nuclear, NMR and resistivity measurements 

PUBLICATION-DATE: July 22, 2004 



INVENTOR- INFORMATION: 
NAME 

Schoen, Juergen S. 
Fanini, Otto N. 
Georgi, Daniel 



CITY 
Leoben 
Houston 
Houston 



STATE 

TX 

TX 



COUNTRY 

AT 

US 

US 



RULE-47 
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Fuf! ] Till == } Ciiition ] Fir*ht Ctest'ti'sZik-rT] t-ate f fM^rK-T] Sec;uersO; ; j /^tfcChiT.eritt: j Cliim-j n!HC 1 I'izw 



□ 3. Document ID: US 20030105590 Al 

L13: Entry 3 of 10 File: PGPB Jun 5, 2003 

PGPUB-DOCUMENT-NUMBER: 2 0030105590 
PGPUB- FILING-TYPE : new 

DOCUMENT- IDENTIFIER : US 20030105590 Al 

TITLE: WATER SATURATION AND SAND FRACTION DETERMINATION FROM BOREHOLE RESISTIVITY 
IMAGING TOOL, TRANSVERSE INDUCTION LOGGING AND A TENSORIAL WATER SATURATION MODEL 

PUBLICATION-DATE: June 5, 2003 



INVENTOR- IN FORMAT I ON : 
NAME 

Mollison, Richard A. 
Schoen, Juergen H. 
Fanini, Otto N. 
Kriegshauser, Berthold 
Pavlovic, Milomir 



CITY 


STATE 


COUNTRY 


Tomball 


TX 


US 


Leoben 


TX 


AT 


Houston 


TX 


US 


Houston 


TX 


us 


Houston 




us 



RULE-47 



US-CL-CURRENT: 702/7 



Titts | f. a;n h j Front J F . ' I l l3^m m.i | I jU \ Raf*f» n< 



□ 4. Document ID: US 20020101235 Al 

L13: Entry 4 of 10 File: PGPB Aug 1, 2002 

PG PUB -DOCUMENT -NUMBER : 20020101235 
PGPUB- FI LING-TYPE : new 

DOCUMENT-IDENTIFIER: US 20020101235 Al 

TITLE: Combined characterization and inversion of reservoir parameters from 
nuclear, NMR and resistivity measurements 

PUBLICATION-DATE: August 1, 2002 



INVENTOR-INFORMATION : 
NAME 

Schoen, Juergen S. 
Fanini, Otto N. 
Georgi, Daniel 



CITY 
Leoben 
Houston 
Houston 



STATE 

TX 

TX 



COUNTRY 

AT 

US 

US 



RULE-47 



US-CL-CURRENT: 32 4/303 
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□ 5. Document ED: US 6899137 B2 

L13: Entry 5 of 10 



File: USPT 



May 31, 2005 



US- PAT-NO: 68 99137 

DOCUMENT-IDENTIFIER : US 6899137 B2 

TITLE: Microf abricated elastomeric valve and pump systems 
DATE-ISSUED: May 31, 2005 



INVENTOR- IN FORMAT I ON : 
NAME 

Unger; Marc A. 
Chou; Hou-Pu 
Thorsen; Todd A. 
Scherer; Axel 
Quake; Stephen R. 
Enzelberger; Markus M. 
Adams; Mark L. 
Hansen; Carl L. 



CITY 

South San Francisco 

Foster City 

Pasadena 

Laguna Beach 

San Marino 

Pasadena 

Pasadena 

Pasadena 



STATE ZIP CODE COUNTRY 

CA 

CA 

CA 

CA 

CA 

CA 

CA 

CA 



US-CL-CURRENT: 137/833; 137 / 565. 16 , 137/803, 137/807, 137/824 



□ 6. Document ID: US 671 1502 B2 

L13: Entry 6 of 10 



File: USPT 



Mar 23, 2004 



US-PAT-NO: 6711502 

DOCUMENT-IDENTIFIER: US 6711502 B2 

TITLE: Water saturation and sand fraction determination from borehole resistivity 
imaging tool, transverse induction logging and a tensorial water saturation model 

DATE-ISSUED: March 23, 2004 



INVENTOR-INFORMATION: 
NAME 

Mollison; Richard A. 
Schoen; Juergen S. 
Fanini; Otto N. 
Kriegshauser; Berthold 
Pavlovic; Milomir 



CITY 

Tomball 

Leoben 

Houston 

Houston 

Houston 



STATE 
TX 

TX 
TX 
TX 



ZIP CODE 



COUNTRY 



AU 



US-CL-CURRENT: 702/6; 702/2, 702 /7 
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□ 7. Document ID: US 6686736 B2 
L13: Entry 7 of 10 



File: USPT 



Feb 3, 2004 



US- PAT-NO: 668 673 6 

DOCUMENT-IDENTIFIER: US 6686736 B2 

TITLE: Combined characterization and inversion of reservoir parameters from 
nuclear, NMR and resistivity measurements 

DATE-ISSUED: February 3, 2004 



INVENTOR- INFORMATION : 
NAME 

Schoen; Juergen S. 
Fanini; Otto N. 
Georgi; Daniel 



CITY 
Leoben 
Houston 
Houston 



STATE 

TX 
TX 



ZIP CODE 



COUNTRY 
AT 



US-CL-CURRENT : 324/303 



□ 8. DocumentID: US 6493632B1 

L13: Entry 8 of 10 



File: USPT 



Dec 10, 2002 



US-PAT-NO: 6493632 

DOCUMENT-IDENTIFIER: US 6493632 Bl 

TITLE: Water saturation and sand fraction determination from borehole resistivity 
imaging tool, transverse induction logging and a tensorial water saturation model 

DATE-ISSUED: December 10, 2002 



INVENTOR-INFORMATION: 
NAME 

Mollison; Richard A. 
Schoen; Juergen S. 
Fanini; Otto N. 
Kriegshauser; Berthold 
Pavlovic; Milomir 



CITY 

Tomball 

Leoben 

Houston 

Houston 

Houston 



STATE 
TX 

TX 
TX 
TX 



ZIP CODE 



COUNTRY 



AT 



US-CL-CURRENT: 702/2; 703/10 



Citation . Front 
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□ 9. Document ED: US 6470274 Bl 

L13: Entry 9 of 10 File: USPT Oct 22, 2002 

US-PAT-NO: 6470274 

DOCUMENT-IDENTIFIER: US 6470274 Bl 

TITLE: Water saturation and sand fraction determination from borehole resistivity 
imaging tool, transverse induction logging and a tensorial dual water saturation 
model 



DATE-ISSUED: October 22, 2002 



INVENTOR- INFORMATION: 
NAME 

Mollison; Richard A. 
Fanini; Otto N. 
Kriegshauser; Berthold 
Pavlovic; Milomir 



CITY 

Tomball 

Houston 

Houston 

Houston 



STATE 

TX 

TX 

TX 

TX 



ZIP CODE 



COUNTRY 



US-CL-CURRENT: 702/7; 702/12 



;<nt | Review if Cl3f~frio3tivf:~| Vsl* f R*f*fcfivi 



File: USPT 



Jul 3, 2001 



□ 10. Document ID: US 6255819 Bl 

L13: Entry 10 of 10 
US- PAT-NO: 6255819 

DOCUMENT-IDENTIFIER: US 6255819 Bl 



TITLE: System and method for geologically-enhanced magnetic resonance imaging logs 
DATE-ISSUED: July 3, 2001 



INVENTOR-INFORMATION: 
NAME 

Day; Peter Ian 
Tutunji; Tarek A. 
Hagiwara; Teruhiko 

US-CL-CURRENT: 324/303; 324/300 



CITY STATE ZIP CODE 

Houston TX 

Houston TX 

Houston TX 



COUNTRY 



■C lass Ma Worn &; £• a te -i ■ R *f * r c 



Generate Collection 




| Print I 


Fwd Refs j 


Bkwd Refs 





Terra 



MAGNETIC 





Generate 


lllll 





Documents 



1396880 
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MAGNETICS 


13310 


RESONANCE 


284416 


RESONANCES 


16359 


MRI 


28372 


MRIS 


369 


NMR 


149933 


NMRS 


236 


(12 AND (MRI OR (MAGNETIC ADJ RESONANCE) OR 
NMR)).USPT ) PGPB,JPAB > EPAB,DWPI,TDBD. 


10 


(LI 2 AND ((MAGNETIC ADJ RESONANCE) OR MRI OR 
NMR)).USPT,PGPB,JPAB,EPAB,DWPI,TDBD. 


10 




Display Format: 1- : l !* an 9 e For,T,at :: j 



Previous Page Next Page Go to Doc# 
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Generate 


Collection 






Print 


j 




Fwd Refs I 


Bkwd Refs 






Generate OACI 





Search Results - Record(s) 1 through 5 of 5 returned. 



□ 1. Document ID: US 4661234 A 
Using default format because multiple data bases are involved. 



L42: Entry 1 of 5 



File: USPT 



Apr 28, 1987 



US-PAT-NO: 4661234 

DOCUMENT- IDENTIFIER : US 4661234 A 

TITLE: Air-fuel ratio sensor and apparatus using the same 
DATE-ISSUED: April 28, 1987 



INVENTOR- INFORMATION: 
NAME 

Takahashi; Hideaki 
Kondo; Haruyoshi 
Takeuchi; Takashi 
Hayakawa; Kiyoharu 



CITY 

Aichi 

Aichi 

Aichi 

Aichi 



STATE ZIP CODE 



COUNTRY 

JP 

JP 

JP 

JP 



US-CL-CURRENT: 204/406; 204/412, 204/425, 204/426, 338/34 



Classification 



□ 2. Document ID: US 3834943 A 

L42: Entry 2 of 5 File: USOC 

US-PAT-NO: 3834943 

DOCUMENT-IDENTIFIER: US 3834943 A 



Sep 10, 1974 



TITLE: ELECTROLYTE-ELECTRODE UNIT FOR SOLID-ELECTROLYTE FUEL CELL AND PROCESS FOR 
THE MANUFACTURE THEREOF 

DATE-ISSUED: September 10, 1974 

INVENTOR-NAME: TANNENBERGER H; VAN DEN BERGHE P 
US-CL-CURRENT: 429 /33, 429/41 











: Review j 
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□ 3. Document ID: US 3479581 A 

L42: Entry 3 of 5 File: USOC 

US- PAT-NO : 3479581 

DOCUMENT-IDENTIFIER: US 3479581 A 



Nov 18, 1969 



TITLE: VERTICAL RESISTIVITY LOGGING BY MEASURING THE ELECTRIC FIELD CREATED BY A 
TIME-VARYING MAGNETIC FIELD 

DATE-ISSUED: November 18, 1969 

INVENTOR-NAME: RUNGE RICHARD J 

US-CL-CURRENT: 324/338; 324/366, 324 / 367 



Foil ]~ Title f Citation [ riord j Review ] edification [ Date ] 



May 29, 1962 



□ 4. Document ID: US 3037105 A 

L42: Entry 4 of 5 File: USOC 

US-PAT-NO: 3037105 

DOCUMENT- IDENTIFIER : US 3037105 A 

TITLE: Methods and apparatus for the induction welding of tubing 
DATE-ISSUED: May 29, 1962 
INVENTOR-NAME: FRED KOHLER 
US-CL-CURRENT: 219/612; 219/613 



i | C rtativn | ; Front |: Riyisw i C \a< s ftfratfan : |;; D ate : 



□ 5. Document ID: US 2784349 A 

L42: Entry 5 of 5 File: USOC 

US-PAT-NO: 2784349 

DOCUMENT-IDENTIFIER: US 2784349 A 
TITLE: Electric arc welding 
DATE-ISSUED: March 5, 1957 
INVENTOR-NAME: ANDERSON NELSON E 

US-CL-CURRENT: 315/176; 123/606, 219 / 130.4 , 315/257 



Mar 5, 1957 
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\? \ :G\ear ,\ Generate Collection J Print \ Fwd Refs. | Bkwd Refs j Generate OACS \ 



Term 


Documents 


PARALLEL 


3307749 


PARALLELS 


17826 


PERMEABILITY 


216475 


PERMEABILITffiS 


8660 


PERMEABILITYS 


4 


RATIO 


2313132 


RATIOS 


438098 


H0RIZ0NTALS3 


0 


HORIZONTAL 


2068539 


HORIZONTALA 


32 


HORIZONTAL AD 


1 


(L37 AND ((H0RIZ0NTAL$3 OR VERTICALS3 OR 
TRANSVERSS5 OR L0NGITUDS4 OR PARALLEL OR 
PERPENDICULARS2 OR LONGITUDINAL$2) WITH 
(RESISTIVS5) WITH (PERMEABILITY OR PERMEABLS3) 
WITH(RATIO))).PGPB,USPT,USOC,EPAB,JPAB,DWPI,TDBD. 


5 



There are more results than shown above. Click here to view the entire set. 



y.v.v.v.v.v.v.v.v 



Display Format : j- 
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Hit List 





Clear 






Generate 


Collection 


n 




Print 1 




Fwd Refs 










Generate OACS 


i 





Search Results - Record(s) 1 through 6 of 6 returned. 



□ 1. DocumentED: US 20040140801 Al 
Using default format because multiple data bases are involved. 

L47: Entry 1 of 6 File: PGPB 



Jul 22, 2004 



PGPUB-DOCUMENT-NUMBER: 2 00401408 01 
PGPUB- FILING-TYPE : new 

DOCUMENT-IDENTIFIER: US 20040140801 Al 

TITLE: Combined characterization and inversion of reservoir parameters from 
nuclear, NMR and resistivity measurements 

PUBLICATION- DATE : July 22, 2004 



INVENTOR-INFORMATION: 
NAME 
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DOCUMENT-IDENTIFIER: US 5463549 A 

TITLE: Method and apparatus for determining permeability of subsurface formations 
Abstract Text (1) : 

Characteristics of invasion profiles exhibiting effects of gravity, and 
particularly their manifestations on resistivity logs of a plurality of different 
radial depths of investigation, are used to determine formation characteristics, 
including vertical permeability . One embodiment includes the following steps: 
suspending a logging tool in the borehole ; producing a plurality of resistivity 
measurements, having respectively different radial depths of investigation, as the 
logging device is moved through the borehole, to obtain a plurality of resistivity 
logs; determining the presence of a buoyancy marker in the resistivity logs in a 
formation bed invaded with filtrate from the drilling mud, and the extent of the 
buoyancy marker; and determining permeability of the formations as a function of 
the extent of the buoyancy marker. 

Assistant Examiner (1) : 
Thomas ; Joseph 

Brief Summary Text (2) : 

This invention relates to well logging of earth boreholes and, more particularly, 
to a method and apparatus for determining permeability of earth formations 
surrounding a borehole . 

Brief Summary Text (3) : 

It is well known that during the drilling of an earth borehole, such as an oil or 
gas well, fluid in the drilling mud begins invading the surrounding formations, due 
to the difference in pressure between the mud in the well bore and the formation 
fluids. The extent of invasion depends upon various factors, including the nature 
of the mud cake formed on the well bore, and the porosity and permeability of the 
surrounding formations . Invasion by the mud filtrate in a given region of the well 
bore will typically continue for a period of time after the region has been 
drilled. Characteristics of the "invaded zone" are important in the determination 
of formation characteristics, and various types of well logging equipment can 
provide measurements of invaded zone parameters that are useful in formation 
evaluation. 

Brief Summary Text (4) : 

Typically, the invasion profile around the well bore is considered to be generally 
cylindrical in shape. For example, in a particular region (or bed) of the 
formation, the invaded zone is normally visualized as a cylindrical annulus having 
a thickness that can vary from a fraction of an inch to several feet. 

Brief Summary Text (5) : 

Although conventional well logging models assume a cylindrical invasion profile, it 
was recognized decades ago that the actual profile of invading fluid in a formation 
bed can be substantially affected by gravity. H. G. Doll noted that in salt -water 
bearing sands of high permeability, it was often observed that the depth of 
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invasion by mud filtrate is quite small near the bottom boundaries of the beds, 
whereas the top parts of the beds are deeply penetrated. (H. G. Doll, "Filtrate 
Invasion In Highly Permeable Sands ", H. G. Doll, The Petroleum Engineer, 1955) . He 
indicated that this observation was brought to light by resistivity logs, in which 
the shallow investigation log gives the same very low readings as the deep 
investigation log over the lower part of the beds, and an appreciably higher 
reading over the upper part. This observation suggested to Doll that, in this kind 
of formation, after the mud filtrate has been forced through the mud cake, it does 
not flow horizontally into the formation, but has a tendency to rise, since it is 
often less saline (and therefore less dense) than the interstitial water originally 
located in the pores of the invaded formations . As a consequence of the difference 
in densities, the mud filtrate is subject to an ascensional force, and its velocity 
is the resultant of two components: one vertical (upward) component, caused by the 
ascensional force, and one horizontal (radial) component, due to the rate of 
filtration. Doll described the nature of the invasion front as it advances, with 
time, into the formation, and set forth a number of relationships representing the 
radial and vertical movement and positions of the advancing filtrate. 

Brief Summary Text (6) : 

Although Doll's work laid a basic foundation for understanding vertical movement of 
the invading filtrate, succeeding years have produced little, if any, development 
on this basic foundation that has resulted in new and commercially useful 
techniques in the logging art for evaluating earth formations . 

Brief Summary Text (7) : 

Determination of formation permeability, including the vertical and horizontal 
components thereof, is important in evaluating the behavior. and producibility of a 
formation . However, it is generally considered difficult to obtain accurate 
permeability information, particularly from conventional types of logging 
measurements. 

Brief Summary Text (8) : 

It is among the objects of the present invention to provide improved logging 
techniques and apparatus for evaluating earth formations surrounding a borehole to 
determine characteristics that are affected by the gravity-induced movement and 
location of invading fluids,. It is also among the objects of the invention to 
provide improved technique and apparatus for determining permeability of earth 
formations . 

Brief Summary Text (10) : 

In accordance with an aspect of the present invention, characteristics of invasion 
profiles exhibiting effects of gravity, and particularly their manifestations on 
resistivity logs of a plurality of different radial depths of investigation, are 
used to determine formation vertical permeability . 

Brief Summary Text (11) : 

An embodiment of a method in accordance with an embodiment of the invention for 
determining .the vertical permeability of formations surrounding a- borehole drilled 
in the earth using drilling fluid, includes the following steps: suspending a 
logging tool in the borehole ; producing a plurality of resistivity measurements 
(preferably at least three) , having respectively different radial depths of 
investigation, as the logging device is moved through the borehole, to obtain a 
plurality of resistivity logs; determining the presence of a buoyancy marker in the 
resistivity logs in a formation bed invaded with filtrate from the drilling mud, 
and the extent of the buoyancy marker; and determining permeability of the 
formations as a function of the extent of the buoyancy marker. In a preferred 
embodiment of the invention, there is derived the elapsed time between the drilling 
of the portion of the borehole for which permeability is being determined and the 
logging of said resistivity measurements in said portion of the borehole, and 
permeability is determined as a function also of said elapsed time. In this 
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embodiment, the buoyancy marker results from a curved segment in the invasion 
profile of the bed, and the extent of the marker is determined from the extent of a 
horn in the resistivity log which has the relatively shallowest radial depth of 
investigation of said resistivity logs. Also in this embodiment, the logging tool 
further includes means for generating a log indicative of bed boundaries, and the 
step of determining a buoyancy marker in the resistivity logs comprises locating 
said horn at one of the bed boundaries. 

Brief Summary Text (12) : 

In accordance with a further feature of the invention, the acceptability of 
provisionally determined values of formation permeability, such as the vertical 
component thereof, is determined by generating a model invasion profile, generating 
simulated resistivity logs using the model invasion profile, and accepting the 
provisional value of permeability if the simulated logs compare favorably with the 
measured logs. 

Drawing Description Text (3) : 

FIGS. 2A-2D illustrate the development of an invading filtrate front as a filtrate 
invades a bed in a situation where the invading filtrate is less dense than the 
formation fluid. 

Drawing Description Text (4) : 

FIG. 3 illustrates an invasion profile in a situation where the invading filtrate 
is denser than the formation fluid, and before steady state is reached. 

Drawing Description Text (5) : 

FIG. 4 illustrates an invasion profile in a situation where the invading filtrate 
is denser than the formation fluid, after steady state is reached. 

Drawing Description Text (7) : 

FIG. 5B illustrates three resistivity logs, at different radial depths of 
investigation, of a type that would be obtained when logging the FIG. 5A formation . 



Drawing Description Text (8) : 

FIG. 5C illustrates several resistivity logs, at different radial depths of 
investigation, of a type that would be obtained when logging the FIG. 5A formation . 

Drawing Description Text (9) : 

FIG. 5D illustrates representative static potential' and gamma ray logs of the type 
that would be obtained when logging the FIG. 5A formation . 

Drawing Description Text (10) : 

FIG. 6A illustrates an invasion profile in a permeable bed with a very thin shale 
streak, invaded by filtrate that is less dense than the formation fluid. 

■ Drawing Description Text (11): . . . 

FIG. 6B illustrates three resistivity logs, having different radial depths of 
investigation, of the type that would be obtained when logging the formation of 
FIG. 6A. ' 

Detailed Description Text (2) : 

FIG. 1 shows an apparatus which can be used to practice an embodiment of the 
invention for investigation earth formations 31 traversed by a borehole 32, which 
is filled with a drilling fluid 33. The investigating apparatus or logging device 
40 is suspended in the borehole 32 on an armored cable 35, the length of which 
substantially determines the depth of the downhole device 40. The cable length is 
controlled by conventional means at the surface (not shown) . Circuitry 51, shown at 
the surface, although portions thereof may be downhole, represents processing 
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circuitry for the various logging units of apparatus 40. A sheave wheel 96 can be 
conventionally provided for determination of depth of the logging apparatus, and 
signals therefrom are coupled with a recorder 90 which generically represents known 
graphical, electrical, and/or other storage and recording functions performed on 
signals received from processing circuitry 51 and from a computing module 100. The 
computing module 100, which in the illustrated embodiment receives inputs from the 
processing circuitry 51, can be implemented, for example, by a suitably programmed 
general purpose computer, such as a VAX, sold by Digital Equipment Corp. It will be 
understood, however, that a suitable special purpose digital or analog computer, 
which performs functions as described herein, could alternatively be employed. ' 

Detailed Description Text (3) : 

The downhole device 40 comprises a tool string that includes several tools. In the 
present embodiment, the reference numeral 41 represents one or more resistivity 
tools which collectively include the capability of measuring resistivity at a 
plurality of radial depths of investigation and, preferably, at least three radial 
depths of investigation. It is well known in the art that one or more resistivity 
logging devices can be employed for this purpose. [Resistivity and conductivity are 
inverses, and it will be understood that references herein to resistivity imply 
that conductivity could be substituted, in an inverse sense.] A single device which 
obtains measurements at several radial depths of investigation is the so-called 
Array Induction Tool or "AIT", which employs an array of induction coils and signal 
processing in the obtainment of its resistivity measurements. Reference can be 
made, for example, to the U.S. Pat. No. 4,873,488 and to Hunka et al . , "A New 
Resistivity Measurement System For Deep Formation Imaging And High Resolution 
Formation Evaluation", SPE-20559, In Annual Technical Conference And Exhibition 
Proceedings, Formation Evaluation and Reservior Geology, Society of Petroleum 
Engineers, 1991. Another technique for obtaining resistivity measurements for at 
least three different depths of investigation is to utilize measurements from a so- 
called "DIT-SFL" (Dual Induction Tool — Spherically Focused Log) apparatus to obtain 
resistivity measurements at a shallow radial depth of investigation ("R.sub.SFL 11 ) 
of about 16 inches, and at intermediate and deep radial depths of investigation 
("R.sub.ILM " and "R.sub.ILD ", respectively) of about 30 inches and 50 inches, 
respectively. [The terms "shallow", "intermediate", and "deep", as used herein, are 
relative terms.] Reference can be made to T. D. Barber, "Introduction To The Phasor 
Dual Induction Tool", Journal Of Petroleum Technology, September 1985; N. Schuster 
et al., "Application of the IFF/Sonic Combination Tool To Gulf Coast Formations ", 
Gulf Coast Assoc. of Geological Societies Transactions, 1971; and Barber, "Phasor 
Induction Tool", Schlumberger WTS, 1989. It will be understood, however, that other 
resistivity logging device (s) can alternatively be employed. 

Detailed Description Text (5) : 

To facilitate understanding of the operation hereof, consider first an example, 
illustrated in FIGS. 2A-2D, of a model formation of the type treated by H. G. Doll 
(supra), wherein a vertical well bore traverses a single homogeneous horizontal 
bed, bounded by impermeable formations (e.g. shales) , and having horizontal 
permeability, k.sub.h, and vertical permeability, k.sub.v. It is assumed that the 
well bore appears instantaneously in the bed at the time of drilling (t=0) , with a 
mudcake in "dynamic equilibrium" with the conditions in the. well, which implies 
that the invading filtrate enters the bed at a constant rate. This is good 
approximation of the generally understood dynamics of mudcake formation and 
invasion. [The volume of filtrate which invades the bed before a mud cake develops, 
is relatively small, and the amount of filtrate which invades the bed after static 
conditions prevail in the well is also relatively small.] 

Detailed Description Text (6) : 

The radial (Darcy) speed at which the filtrate enters the formation, u.sub.f, is 
assumed to be controlled by the properties of the mud. [The volume of filtrate per 
unit time per unit distance entering the formation is u.sub.f times the wellbore 
perimeter.] This is a direct consequence of the relatively large hydrodynamic 
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resistance encountered by filtrate as it flows thorough the mud cake (as compared 
to the formation ) . This presumes the presence of a mud cake, implying k.sub.h of 
the bed is above some minimum value, typically at least 10 mD. Once the horizontal 
permeability exceeds this value (and a mud cake is present) , then the specific 
value of u.sub.f is independent of k.sub.h. 

Detailed Description Text (7) : 

It is also assumed that the invasion process can be modeled by a saturation "step 
profile"; i.e., the filtrate saturation is a constant throughout the invaded zone, 
the annulus of connate water (situated between the filtrate and the formation 
fluid) is not included, and the effect of capillary forces or viscous instabilities 
can be ignored resulting in a transition zone (region between the filtrate and the 
formation fluid) having zero thickness. Rigorously speaking, these conditions model 
aqueous filtrate invading a water zone; however, they may also be appropriate in 
situations when the filtrate and the formation fluid consist of different phases. 
(See e.g. Singer et al . , "The Effect Of Transition Zones On The Response Of 
Induction Logs", SWPA Annual Logging Symposium, 1988, and Allen et al . , Invasion 
Revisited, Oilfield Review, 1991) . 

Detailed Description Text (8 ) : 

FIGS. 2A-2D illustrate the development of the filtrate front (see H. G. Doll, 
supra) as the filtrate invades the bed, in a situation where the invading filtrate 
is less dense than the formation, fluid. FIG. 2A represents time t=0, when invasion 
first starts, and identifies four filtrate material points in the wellbore wall as 
a, b, c, and d. FIGS. 2B, 2C, and 2D show the filtrate front, and the progress of 
the points a, b, c and d, at successively greater times, t.sub.l, t.sub.2 and 
t.sub.3, respectively. Depending on the parameters involved, t.sub.l, t.sub.2 and 
t.sub.3 may be, for example, about one day, five days, and ten days, respectively, 
after drilling. By the time t.sub.3, the filtrate points b, c, and d have entered 
the upper portion of the bed, with the point a about to enter. Time t.sub.3, in 
this example, is approximately the time, T, treated further below, that it takes 
for the vertical layer of the filtrate front to reach a steady state. 

Detailed Description Text (9) : 

The difference in densities between the filtrate and the formation fluid can have a 
positive or a negative value. [In FIGS. 2A-2D the filtrate is less dense than the 
formation fluid (e.g. fresh water filtrate invading saline formation fluid), so the 
filtrate rises, but if the filtrate were more dense than the formation fluid (e.g. 
water filtrate invading hydrocarbon formation fluid), the filtrate would fall.] 
Thus, the accumulate filtrate near a horizontal impermeable bed boundary will be 
referred to as the " horizontal layer", and the remainder of the filtrate front will 
be referred to as the vertical layer. [Initial reference can be made to FIGS. 3 and 
4, which represent the case where the invading filtrate is denser than the 
formation fluid, and are upside-down approximate renditions of FIGS. 2C and 2D, 
respectively.] It is know (e.g. H. G. Doll, supra) that the velocity of the 
filtrate, u.sub.fil, within the formation is given by ##EQU1## where R.sub.wb 
denotes the wellbore radius, and r and z denote, respectively, local unit vectors 
in the radial and axial., direct ions..- o.f , a cylindrical coordinate .system (r,,,..phi., z) 
whose z-axis coincides with the center of the wellbore and points upward (the plane 
z = 0 indicating the location of the lower impermeable bed boundary) . The constant, 
W, is defined by 

Detailed Description Text (10) : 

where .rho. and .mu. denote, respectively, the density and the viscosity of the 
fluid (the subscripts "fil" and "for" referring to the filtrate and the formation 
fluid, respectively), and g denotes the gravitational constant (980 cm/sec . sup . 2 ) . 
The radial component of filtrate velocity is identical to the case when the 
filtrate has the same density as the formation fluid. Its vertical component is a 
constant, W, not varying either with position within the formation or over time. 
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Detailed Description Text (11) : 

The filtrate velocity of equation (1) can be used to locate the filtrate points of 
FIGS. 2A-2D at the four successive times t.sub.l, t.sub.2, and t.sub.3. At each 
time, all of the points have advanced into the formation the same distance, 
maintaining identical radial positions, and have risen vertically by the same 
amount, retaining their same vertical spacing until they have entered the upper 
portion of the bed. The "curved segment" of the vertical layer of the filtrate 
front, which begins to develop along its lower end, denotes the location of 
filtrate points which have entered the formation at position (r, z) = (R.sub.wb, 0) 
at times later than t=0 (the . phi . -coordinate is not specified due to azimuthal 
symmetry) . The vertical straight-line segment denotes the location of all the 
filtrate points at the wellbore wall at t=0. At time . phi . H . sub . f or /W, denoted by 
T, the shape of the filtrate front has reached its steady state configuration, 
where .phi. denotes the porosity and H. sub. for denotes the bed thickness. This 
corresponds to the approximate time it takes the filtrate point located at 
(R.sub.wb, 0) at time t=0 to travel to the upper part of the bed. 

Detailed Description Text (12) : 

The location of the filtrate front vertical layer, prior to steady state (i.e. 
t .ltoreq. T) , is ##EQU2## where . DELTA.. rho. is the filtrate density minus the 
formation fluid density and, as noted above, R.sub.wb is the radius of the 
wellbore, and u.sub.f is the radial speed of the filtrate as it enters the bed at 
the wellbore wall (as represented by equation (1)), and t is the time of 
observation. The time to steady state, T, is 

Detailed Description Text (13) : 

The equation (3) applies to the "curved segment" of the vertical layer which, in 
FIG. 3, is adjacent to the top of the permeable bed [ ( 1-t . sub . log /T)H.sub.for 
<z<H.sub.for in FIG. 3], and has a vertical extent, H.sub.v. Below the curved 
segment [ 0<z< ( 1-t . sub . log /T)H.sub.for in FIG. 3], the vertical layer is 
represented by ##EQU3## After steady state (t .gtoreq. T) , the equation (3) applies 
to the entire vertical layer (see steady vertical layer illustration of FIG. 4) . 

Detailed Description Text (14) : 

Referring again to the invasion profile prior to steady state (t .ltoreq. T) of 
FIG. 3, the vertical extent of the curved segment of the vertical layer is given by 
##EQU4## ~ ' " " 

Detailed Description Text (15) : 

Applicants have determined that, for t .ltoreq. T, the horizontal and vertical 
extents of the horizontal layer can be represented as: ##EQU5## and for the case 
for t .gtoreq. T, the horizontal and vertical extents of the horizontal layer can 
be represented as: ##EQU6## The polynomial expressions are valid for 0 . 05< . eta . <1, 
and the above solutions were obtained by numerical solution of partial differential 
equations representing the dynamic geometry of the horizontal layer. If desired, 
the location of the filtrate front over the entire permeable zone can be 
constructed by combining the above expressions for the vertical and horizontal 
layers . ...It will , be understood that the. effects described above . are .significant. .-when 
"buoyancy" is relatively large; that is, when the vertically induced speed W is 
greater than the radial filtrate speed at the wellbore wall,, u.sub.r. 

Detailed Description Text (16) : 

In accordance with a feature of the present invention, the presence of a curved 
segment of a vertical layer, in an invasion profile with buoyancy, is detected 
using resistivity logging measurements, and the vertical extent, H.sub.v of the 
curved segment is determined. The vertical permeability of the formation in a bed 
of interest can then be determined as a function of H.sub.v, the porosity .phi., 
.DELTA. . rho. , .mu . . sub . f il, and t.sub.log, which is the time after drilling that 
the resistivity measurement were taken. The vertical permeability can be determined 
from ##EQU7## which follows from equation (6) above. 
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Detailed Description Text (17) : 

In accordance with a further feature of the invention, the acceptability of 
provisionally determined values of formation permeability, such as the vertical 
component thereof, is determined by generating a model invasion profile, generating 
simulated resistivity logs using the model invasion profile, and accepting the 
provisional value of permeability if the simulated logs do not differ substantially 
from the measured logs. 

Detailed Description Text (18) : 

Referring to FIGS . 5A-5D, there are shown examples of representative log traces of 
the types obtained from the resistivity logging device (s) of the FIG. 1 embodiment 
in a permeable formation bed between shale shoulder beds in a case where the 
filtrate is lighter than the fluid in the invaded formations . In this example, 
fresh water filtrate invades a water zone with R.sub.t <R.sub.xo <R. sub . shale . The 
approximate invasion profile is shown in FIG. 5A, and it is seen to have the 
characteristic pattern described in conjunction with FIGS. 2 and 3, with the shaded 
invasion region exhibiting a curved segment 510, a vertical straight line segment 
520, and an accumulated ( horizontal ) layer 530. The horizontal layer is adjacent 
the top of the bed and the curved segment is adjacent the bottom of the bed, which 
follows from the buoyancy of the lighter invading filtrate for this example. FIG. 
5D illustrates the SP and GR logs for the FIG. 5A example, and can be utilized, 
inter alia, for bed boundary determination. FIG. 5B illustrates the log traces, as 
a function of depth, for the DIT-SFL, which produces the shallow (radial depth of 
investigation) resistivity R.sub.SFL, the intermediate resistivity R.sub.ILM, and 
the deep resistivity R.sub.ILD. FIG. 5C shows the AIT logs at radial depths of 
investigation of 10, 20, 30, 60 and 90 inches. In this example, the depth of 
invasion at the vertical straight line segment is about 15 inches at the 
illustrated time. In both the DIT-SFL logs of FIG. 5B and the AIT logs of FIG. 5C, 
the deepest resistivity measurements (R.sub.ILD for FIG. 5B, and R. sub. 60 and 
R. sub. 90 for FIG. 5C) see substantially only the uninvaded formation, and exhibit a 
relatively low resistivity in the permeable water -containing bed. The shallowest 
measurements R.sub.SFL (for FIG. 6) and R. sub. 10 (for FIG. 7) are most strongly 
influenced by the invaded zone in this example, and exhibit resistivities 
characteristic of the fresh water filtrate in the bed. This resistivity is seen to 
be exhibited over the portions of the bed where the vertical straight-line segment 
520 and the accumulated (or horizontal ) layer 530 are located. However, in the 
region of the curved segment 510, there is seen to be a horn in the shallow 
resistivity logs, the horn having a vertical extent which approximately corresponds 
to the vertical extent of the curved segment, H.sub.v. The logs having intermediate 
depths of radial investigation (R.sub.ILM in FIG. 5B, and R. sub. 20 and R. sub. 30 in 
FIG. 5C) , which are influenced to a much smaller extent by the invaded zone, 
exhibit a smaller perturbation in the region of the curved segment. The horizontal 
layer does not have sufficient vertical extent to have a noticeable effect on these 
resistivity logs. However, in conditions where k.sub.h is approximately equal to 
k.sub.v, the horizontal layer may cause a noticeable deflection on shallow 
resistivity logs in the direction of R.sub.xo. Also, if resistivity devices (s) with 
higher vertical resolution are used, the horizontal layer would exhibit ...larger 
effects. In the shale regions, where there is no invasion, all resistivity 
measurements exhibit the same relatively high resistivity . 

Detailed Description Text (19) : 

FIG. 6A illustrates a buoyancy induced invasion profile for a situation where there 
is a very thin shale streak in a permeable bed which, again, is invaded by fresh 
water filtrate. Again, the filtrate is lighter than the formation fluid, and will 
tend to rise, thereby giving rise to a horizontal layer 630 below the shale streak 
and a curved segment 610 above the shale streak. FIG. 6B illustrates the DIT-SFL 
logs for this situation where, again, R.sub.t >R.sub.xo >R . sub . shale . As before, 
the resistivity log for the deepest radial depth of investigation, R.sub.ILD, sees 
the uninvaded formation, and is not affected by the very thin shale steak or the 
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invasion. Also, as before, the shallowest measurement exhibits a horn that results 
from the curved segment (610), most of the vertical extent of the horn being above 
the shale streak. The excursion of the horn goes toward R.sub.t rather than toward 
R. sub. shale The intermediate depth of investigation measurement (R.sub.ILM) 
exhibits a smaller perturbation due to the curved segment 610. Again, it is 
believed that the accumulated ( horizontal ) layer 630 does not noticeabley affect 
these logs. Any significant permeability in the shale would allow the accumulating 
layer of filtrate to pass vertically through it, resulting in deeper invasion in 
the bed above the shale, thus reducing or eliminating the signature in the shallow 
measurement (i.e., the horn). A distinct signature indicates that the radial extent 
of the impermeable shale is at least as large as the accumulated layer. 

Detailed Description Text (20) : 

Referring to FIG. 7, there is shown a flow diagram of a routine for programming the 
processor of computing module 100 (FIG. 1) in accordance with an embodiment of the 
invention. The block 710 represents the storing of logging measurements from the 
tool string of the FIG. 1 embodiment, to obtain a plurality of resistivity and 
other logs, as described. The stored logs can be displayed and/or recorded, in 
conventional fashion, by recorder 90, as represented by the block 715. The block 
720 represents the determination of bed boundaries from the SP and/or GR logs, 
which can be performed in conventional fashion by detecting, for example, the types 
of excursions illustrated in FIG. 5D. 

Detailed Description Text (21) : 

The block 725 represents determination of the presence of a curved segment (of the 
type illustrated, for example, at 510 in FIG. 5A) from resistivity logs having 
relatively shallow and relatively deep radial depths of investigation. In the 
present embodiment, a curved segment is identified by observing a large deflection 
(horn) , at a bed boundary, in the shallow resistivity log toward the resistivity of 
the deep resistivity log. The extent of the horn (in the borehole axial direction), 
from which one can determine H.sub.v, is taken in this embodiment to be the 
distance between the boundary and the position where the horn reaches the bed 
resistivity . The block 730 represents the deriving of information, including 
measurement information and/or information obtained from logging, drilling 
parameters, local geology information, and the like, to obtain (and/or receive, by 
operator input) values for porosity, .phi., the time since logging, t.sub.log, and 
the densities of the fluids involved and the viscosity of the filtrate. For 
example, porosity, as a function of depth, may be obtained from the logging device 
43 of FIG. 1, the fluid properties may be obtained from the known mud properties 
and local geology information, and the time since drilling of the section of 
borehole being considered may be derived by the operator from standard drilling 
record information at the well site. 

Detailed Description Text (22) : 

The block 740 is next entered, this block representing the determination of a 
provisional value of vertical permeability, k.sub.v, in accordance with 
relationship (11) above. The block 745 is then entered, this block representing the 
generation -of an invasion prof ile , of the formations from the^resistiyity logs and--, 
the vertical permeability determined at the block 740. 

Detailed Description Text (23) : 

In the well logging art, so-called "inverse modeling" and so-called "forward 
modeling" are used for obtaining models of formations, or of simulated logging 
device responses, for purposes of evaluating the nature of subsurface formations . 
[Reference can be made to B. Anderson et al . , "ELMOD — Putting Electromagnetic 
Modeling To Work To Improve Resistivity Log Interpretation", Transactions of SPWLA 
30th Annual Logging Symposium, June, 1989; and D. Allen et al . , "Invasion 
Revisited", Oilfield Review, July, 1991; and A. Q. Howard, "A New Invasions Model 
For Resistivity Log Interpretation", The Log Analyst, March-April, 1992.] In 
inverse modeling, one may have measurements from a plurality of logging devices 



http://westbrs:9000/bin/gat^ 9/19/05 



(such as resistivity logs) and, from the logs, develop a model of the formations, 
for example a model which delineates bed and invasion geometries and resistivities . 
Conversely, a forward model for a particular logging tool can begin with a 
formation model and is used to derive a simulated log that would be produced by 
that logging tool traversing a borehole though the formation model. For example, in 
forward modeling of a resistivity tool, such as an induction tool having a number 
of coils, the tool characteristic is readily derived from the physical 
characteristics of the coils, their locations, and the signals applied thereto. In 
the present embodiment, a model, with buoyancy taken into account, is developed 
using equations (3) and (5) , as indicated, and the previously determined 
provisional value of k.sub.v. If desired, a model which does not take buoyancy into 
account (i.e., with a cylindrical invasion profile and no vertical movement of the 
fluids considered), can be developed for comparison purposes. Next, the block 760 
represents generation of simulated logs from the formation model (s) using forward 
modeling. In this example, the simulated logs for the DIT-SFL are developed from 
the previously derived invasion model which takes buoyancy into account. [If 
desired, simulated logs, without buoyancy, can also be generated by forward 
modeling from the model without buoyancy.] The simulated logs (with buoyancy) are 
then compared to the actual logs (decision block 765), and the provisional value of 
k.sub.v is accepted (block 771) if the comparison is favorable, and rejected (block 
772) if it is not favorable. Various criteria can be used for the comparison. For 
example, the presence of the horn in the shallow resistivity logs can be required 
to match within a predetermined percentage in position and amplitude, as a 
condition of acceptance. The next section of the logs can then be processed (block 
780), with the block 710 being re-entered to initialize. 

Detailed Description Text (24) : 

In relatively thin beds, in which the vertical layer has achieved its final shape 
at the time of logging (T<t . sub . log) , the resistivity logs can determine, at most, 
u.sub.f /k.sub.v. This can be seen as follows. Since the vertical layer has 
achieved its final shape, the radial position- of the invasion front is given by 
equation (3). Again, if . mu . . sub . f il, . DELTA. . rho . , R.sub.wb and H. sub. for are 
determined from logs, local knowledge, or laboratory analysis (as above) , then this 
expression for. the position of the front represents a one parameter family of 
shapes, the parameter being the remaining unknown quantity u.sub.f /k.sub.v. The 
value of u.sub.f /k.sub.v can then be determined such as by iteratively forward 
modeling the response of the resistivity tool(s) until the resistivity logs 
generated by the forward modeling calculations approximately matches those on the 
field log. Thus, obtaining the vertical permeability in thin beds requires the 
additional knowledge of u.sub.f, the radial speed of the filtrate as it enters the 
bed. Since u.sub.f depends primarily on the properties of the mud and the shear 
stress on the wellbore wall produced by the mud flowing through the wellbore, and 
not on the formation properties,- then the value of u.sub.f evaluated in thick 
neighboring beds may be used. [In the thick beds, u.sub.f can be determined from 
the resistivity logs associated with the vertical straight-line segment of the 
vertical layer as expressed by equation (5), with the diameter of invasion 
determined such as by inversion modeling the resistivity logs or by iteratively 
t-„f o rwa r d mode Li ng the resistivity, jto ol. (s ) . ] . . . ., w • . j . 

Detailed Description Text (25) : 

Referring to FIG. 8, there is shown a flow diagram of a routine for controiLing the 
processor of computing moduie 100 to obtain a log of vertical permeability values 
using the type of process just described. The portion of the FIG. 7 routine 
described in conjunction with the blocks 710-720 would be the same as previously 
described. The block 830 corresponds to the block 730 of the FIG. 7 routine. The 
block 840 is then entered, this block representing selection of an initial model 
value of u.sub.f /k.sub.v. A conductivity model with an invasion profile that takes 
buoyancy into account is then generated (block 845) , as previously described, using 
a forward model and u.sub.f /k.sub.v in equation (3) to obtain r (radius of the 
invasion front) as a function of z (depth or axial distance) . Simulated logs can 
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then be generated (block 850) using forward modeling, as previously described. The 
decision block 860 is then entered, this block representing the comparison of the 
simulated logs with the actual measured (field) logs. There are various known ways 
in which the comparison can be implemented, such as by computing a -root mean square 
error over a particular region. If the error is above a predetermined threshold 
level, the model value of u.sub.f /k.sub.v is modified (block 865), the block 845 
is re-entered, and the loop 875 continues until the error is less than the 
predetermined error threshold value. The block 880 is then entered, this block 
representing the deriving of u.sub.f, and the storage of the resultant vertical 
permeability, k.sub.v. As described above, the value of u.sub.f can be determined, 
for example, from adjacent beds. The next section of the logs can then be 
considered (block 885), and the block 710 is re-entered for further processing. 

Detailed Description Text (26) : 

The invention has been described with reference to particular preferred 
embodiments, but variations within the spirit and scope of the invention will occur 
to those skilled in the art. For example, while formation beds used for 
illustration have been shown perpendicular to the borehole (zero dip) , it will be 
understood that the invention is equally applicable to dipping formations . 

Other Reference Publication ( 1 ) : 

H. G. Doll, "Filtrate Invasion In Highly Permeable Sands ", The Petroleum Engineer, 
1955. 

Other Reference Publication (2) : 

Hunka et al . , "A New Resistivity Measurement System For Deep Formation Imaging And 
High Resolution Formation Evaluation", SPE-20559, In Annual Technical Conference 
And Exhibition Proceedings, Formation Evaluation and Reservior Geology, Society of 
Petroleum Engineers, 1991. 

Other Reference Publication (4) : 

N. Schuster et al . , "Application of the IFF/Sonic Combination Tool To Gulf Coast 
Formations ", Gulf Coast Assoc. of Geological Societies Transactions, 1971. 

Other Reference Publication (8) : 

B. Anderson et al . , "ELMOD — Putting Electromagnetic Modeling To Work To Improve 
Resistivity Log Interpretation", Transactions of SPWLA 30th Annual Logging 
Symposium, Jun., 1989. 

Other Reference Publication (9) : 

A. Q. Howard, "A New Invasions Model For Resistivity Log Interpretation", The Log 
Analyst, Mar. -Apr., 1992. 

CLAIMS : 

I. A method for determining the permeability of formations surrounding a borehole 
drilled in the earth using drilling fluid, comprising the steps of: 

suspending a logging tool in the borehole ; 

producing a plurality of resistivity measurements, having respectively different 
radial depths of investigation, as said logging device is moved through said 
borehole, to obtain a plurality of resistivity logs; 

determining the presence of a buoyancy marker in said resistivity logs in a 
formation bed invaded with filtrate from said drilling fluid, and the extent of 
said buoyancy marker; and 

determining permeability of said formations as a function of the extent of said 
buoyancy marker. 
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2. The method as defined by claim 1, further comprising deriving of the elapsed 
time between the drilling of the portion of the borehole for which permeability is 
being determined and the logging of said resistivity measurements in said portion 
of the borehole/ and wherein said permeability is determined as a function also of 
said elapsed time. 

3. The method as defined by claim 1, wherein said step of determining permeability 
comprises determining the vertical permeability of said formations . 

4. The method as defined by claim 2, wherein said step of determining permeability 
comprises determining the vertical permeability of said formations . 

5. The method as defined by claim 1, wherein said plurality of resistivity 
measurements, having respectively different depths of investigation, comprises at 
least three resistivity measurements from which at least three logs are obtained. 

6. The method as defined by claim 4, wherein said plurality of resistivity 
measurements, having respectively different depths of investigation, comprises at 
least three resistivity measurements from which at least three logs are obtained. 

7. The method as defined by claim 1, wherein said buoyancy marker comprises a 
curved segment in the invasion profile of said bed, and wherein the extent of said 
marker is determined from the extent of a horn in the resistivity log having the 
relatively shallowest radial depth of investigation of said resistivity logs. 

8. The method as defined by claim 4, wherein said buoyancy marker comprises a 
curved segment in the invasion profile of said bed, and wherein the extent of said 
marker is determined from the extent of a horn in the resistivity log having the 
relatively shallowest radial depth of investigation of said resistivity logs. 

9. The method as defined by claim 6, wherein said buoyancy marker comprises a 
curved segment in the invasion profile of said bed, and wherein the extent of said 
marker is determined from the extent of a horn in the resistivity log having the 
relatively shallowest radial depth of investigation of said resistivity logs. 

10. The method as defined by claim 8, wherein said logging tool further includes 
means for generating a log indicative of bed boundaries, and wherein said step of 
determining a buoyancy marker in said resistivity logs comprises locating said horn 
at one of said bed boundaries. 

11. The method as defined by claim 9, wherein said logging tool further includes 
means for generating a log indicative of bed boundaries, and wherein said step of 
determining a buoyancy marker in said resistivity logs comprises locating said horn 
at one of said bed boundaries. 

12. The method as defined by claim 6, further comprising generating an output log 
of said, determined permeability ^ . .... ....... 

13. A method for determining the vertical permeability of formations surrounding a 
borehole drilled in the earth using drilling fluid, comprising the steps of: 

suspending a logging tool in the borehole ; 

producing a plurality of resistivity measurements, having respectively different 
radial depths of investigation, as said logging device is moved through said 
borehole, to obtain a plurality of resistivity logs; 

determining the presence of a buoyancy marker in said resistivity logs in a 
formation bed invaded with filtrate from said drilling fluid, and the extent of 
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said buoyancy marker; 



determining permeability of said formations as a function of the extent of said 
buoyancy max ke re- 
determining an invasion profile for said bed from said resistivity logs and said 
provisional vertical permeability ; 

generating a plurality of simulated resistivity logs from said invasion profile; 

comparing said simulated resistivity logs with the resistivity logs obtained from 
measurements to determine whether the provisional vertical permeability is 
acceptable; and 

outputting the acceptable vertical permeability value. 

14. The method as defined by claim 13, further comprising deriving of the elapsed 
time between the drilling of the portion of the borehole for which permeability is 
being determined and the logging of said resistivity measurements in said portion 
of the borehole, and wherein said permeability is determined as a function also of 
said elapsed time. 

15. The method as defined by claim 14, wherein said plurality of resistivity 
measurements, having respectively different depths of investigation, comprises at 
least three resistivity measurements from which at least three logs are obtained. 

16. The method as defined by claim 14, wherein said buoyancy marker comprises a 
curved segment in the invasion profile of said bed, and wherein the extent of said 
marker is determined from the extent of a horn in the resistivity log having the 
relatively shallowest radial depth of investigation of said resistivity logs. 

17. The method as defined by claim 15, wherein said buoyancy marker comprises a 
curved segment in the invasion profile of said bed, and wherein the extent of said 
marker is determined from the extent of a horn in the resistivity log having the 
relatively shallowest radial depth of investigation of said resistivity logs. 

18. The method as defined by claim 14, wherein said logging tool further includes 
means for generating a log indicative of bed boundaries, and wherein said step of 
determining a buoyancy marker in said resistivity logs comprises locating said horn 
at one of said bed boundaries. 

19. The method as defined by claim 17, wherein said logging tool further includes 
means for generating a log indicative of bed boundaries, and wherein said step of 
determining a buoyancy marker in said resistivity logs comprises locating said horn 
at one of said bed boundaries. 

20. The method as defined by claim 14, further comprising generating an output log 
of said determined permeability .. , .... .... _ 5 . t . . 

21. Apparatus for determining the permeability of formations surrounding a borehole 
drilled in the earth using drilling fluid, comprising: 

a logging tool suspended in the borehole ; 

means for producing a plurality of resistivity measurements, having respectively 
different radial depths of investigation, as said logging device is moved through 
said borehole, to obtain a plurality of resistivity logs; 

means for determining the presence of a buoyancy marker in said resistivity logs in 
a formation bed invaded with filtrate from said drilling fluid, and the extent of 
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said buoyancy marker; and 

means for determining permeability of said formations as a function of the extent 
of said buoyancy marker. 

22. Apparatus as defined by claim 21, wherein said permeability is determined as a 
function also of the elapsed time between the drilling of the portion of the 
borehole for which permeability is being determined and the logging of said 
resistivity measurements in said portion of the borehole . 

23. Apparatus as defined by claim 22, wherein said logging tool further includes 
means for generating a log indicative of bed boundaries, and wherein said means for 
determining a buoyancy marker in said resistivity logs comprises means for locating 
said horn at one of said bed boundaries. 

24. Apparatus as defined by claim 22, further comprising means for generating an 
output log of said determined permeability . 

25. A method for determining the permeability of formations surrounding a borehole 
drilled in the earth using drilling fluid, comprising the steps of: 

a) suspending a logging tool in the borehole ; 

b) producing at least three resistivity measurements, having respectively different 
radial depths of investigation, as said logging device is moved through said 
borehole, to obtain at least three resistivity logs; 

c) selecting a model value of vertical permeability for a formation bed invaded 
with filtrate from said drilling fluid; 

d) generating from said resistivity logs and said value of vertical permeability, a 
formation resistivity model having an invasion profile which includes the buoyancy 
effect of gravity and depends on said model value of vertical permeability ; 

e) generating simulated logs from said formation resistivity model; 

f ) comparing the simulated logs with the logs obtained from resistivity 
measurements to obtain an error based on dissimilarity between the compared logs; 

g) modifying the model value of vertical permeability ; 

h) repeating steps (d) through (g) until a predetermined criterion of said error is 
met; and 

i) outputting said model value of vertical permeability as the vertical 
permeability in said formation bed. 

26. The method as. defined by. claim 25, wherein said model value ..of vertical ..... 
permeability is in the form of a model ratio of the radial speed of filtrate in 
said bed to vertical permeability in said bed. 

27. The method as defined by claim 26, further comprising repeating the recited 
method steps for other beds, and generating a log of vertical permeability . 
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